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, Abstract. Chirality, i.e., the right- and left-handedness of structure, is one of the key concepts 

in many fields of science including biology, chemistry and physics, and its manipulation is 
an issue of vital importance. The electron spins in solids can form chiral configurations. In 
^ , perovskite manganites i^MnOs {R=Th, Dy,...etc), the Mn-spins form a cycloidal structure, 

which induces ferroelectric polarization (P) through the relativist ic spin-orbit interaction. This 

S magnetism-induced ferroelectricity (multiferroics) and associated infrared-active spin waves 
(electromagnons) open a promising route to control the spins by purely electric means in a 
' very short time. In this paper, we show theoretically with an accurate spin Hamiltonian for 

, TbMnOs that a picosecond optical pulse can switch the spin chirality by intensely exciting the 

Q ■ electromagnons with a terahertz frequency. 

> 

^ ■ 1. Introduction 

I Some structures are not invariant upon the mirror operation just like the right and left hands, 

and are called chiral. This concept is of fundamental importance in chemistry and biology. The 
selective chemical synthesis of chirality (called asymmetric synthesis) is crucial for reliable and 
safe productions of medicines and foods. Note that the chirality of molecules and crystals is 
hard to change once synthesized because it is determined by the chemical bonds. In physics, 
on the other hand, spins in magnets sometimes form a chiral order. The spin chirality is much 
J> I more flexible since the spin is a quantum-mechanical and relativist ic object, which offers an 

opportunity to control the chirality by external parameters. 

The representative laboratory for this spin chirality is the multiferroic transition-metal 
oxides [1]. A model taking into account the spin-orbit interaction predicts the electric 
polarization pij produced by the chiral spin configuration [21 [3] as 

Pij oc Cij X {Si X Sj), (1) 

with Cij being the vector connecting zth and jth sites. Here the vector product Si x Sj is 
the spin chirality, which characterizes a direction of the spin rotation as one goes from i to j. 
This model explains the multiferroic behavior in i^MnOa where the cycloidal Mn-spin order 
induces the ferroelectric polarization P [4J. For example, in TbMnOa [see Figs. [l](a) and (b)], 
the clockwise 6c-plane spin spiral propagating along the +6 direction induces P||+c, while the 
counterclockwise one induces c OlSj. For the a, b and c axes, we adopt the Pbnm setting. 
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Figure 1. (a) Spin configuration and ferroelectric polarization P in the clockwise 6c-plane 
spin spiral and (b) those in the counterclockwise one. (c) Crystal structure, spin exchanges, 
Dzyaloshinskii-Moriya vectors d^j, and local axes attached to the zth MnOe octahedron in 
i?Mn03. Here FM (AFM) denotes (anti) ferromagnetic exchange, (d) Modulation of the in- 
plane ferromagnetic exchanges under £^||±a, which induces slight spin rotations indicated by 
red (blue) arrows under £^||+a (E\\—a). Upper (lower) signs in front of AJab (>0) correspond 
to the modulations under E\\+a {E\\—a). (e) Calculated electromagnon optical spectrum. Inset 
shows the experimental spectrum for TbMnOa from Ref. [13j. 

When the electric polarization is driven by the spin order, it is naturally expected 
that collective excitations of spins (magnons) have an optical activity [3 El Ej. Indeed, 
strong optical absorptions were experimentally observed in i?Mn03 at terahertz frequencies, 
and they were ascribed to magnons activated by the electric- field component of light i.e., 
electromagnons [SI El [IDl [HI [El [13 [HI [E]- In the early stage, the corresponding magnon 
modes were interpreted as rotation of the spin-spiral plane accompanied by oscillation of pij 
in Eq. (1) [0]. However, this interpretation contradicts to the experimental observation about 
the selection rule in terms of the light polarization [TT1[T2]. Afterwards, it turned out that the 
conventional magnetostriction mechanism where pij is given by 

Pij = 7rij(Si • Sj) (2) 

is relevant to the electromagnons in i^MnOs [161 [TT]. Here the vector TVij is nonzero because of 
the orthorhombic lattice distortion without inversion symmetry at the center of the Mn-O-Mn 
bond. The puzzling electromagnon optical spectrum with two specific peaks was successfully 
explained by this mechanism [T7) . 

Under this circumstance, the photo-induced phenomena become a challenging issue for the 
next step. Since pij in Eq. (2) does not require the spin-orbit interaction, its magnitude is much 
larger than that of pij in Eq. (1), which enables us to realize the intense optical excitations 
of magnons. In fact, the real-time dynamics of the photo-excitation is usually a notoriously 
complicated problem since the electronic excitations are the genuine many-body problem of 
various degrees of freedom. Therefore, a reliable theoretical analysis has been almost impossible 
thus far. However, we are now ready to attack such phenomena in i^MnOa theoretically. 
This is because we know that the optical pulse activates mostly the spin systems only via 
the magnetostriction given in Eq. (2), which allows us to neglect electronic excitations at 
higher energies (>1.5 eV). In addition, we have an accurate spin Hamiltonian, which describes 
competitions among various phases in i^MnOs [18], so that the optical switching among them 
can be simulated in a reliable way. Recent rapid experimental progress in intense THz-pulse 



generation also gives a growing importance to the theoretical study on this subject [T9t [201 [2T]. 



2. Model and Method 

To describe the Mn-spin system in TbMnOa, we employ a classical Heisenberg model on a cubic 
lattice, which is given by 



<i,j> i i 

ab 

<hj> <hj> 

Here ix^ iy, iz represent the integer coordinates of the ith Mn ion with respect to the x, y and z 
axes. The first term describes the spin-exchange interactions as shown in Fig. [IJc). The second 
and the third terms stand for the single-ion anisotropy defined with local axes, 77^ and 
attached to each MnOe octahedron tilted in the orthorhombic lattice [18j. The fourth term 
denotes the DM interaction where the DM vectors dij are expressed using five DM parameters, 
c^abi f^ab^ lab^ ^ud j^c^ as given in Ref. ^22j. The last term represents the biquadratic interaction 
originating from the spin-phonon coupling [23]. The Monte-Carlo analysis of this model has 
successfully reproduced the rich phase diagram of i^MnOs [IH]. We perform calculations using 
following parameters: J^5=-0.74, J5=0.64, Jc=1.0, (Q^a6,/3a6,7a6)=(0.1, 0.1, 0.14), (ac,^c)=(0.48, 
0.1), i?=0.2, £"=0.25, and 5biq=0.025, where the energy unit is meV. This parameter set gives 
the 6c-plane transverse spin spiral propagating along the h axis with a wave number g^^O.Svr at 
low temperatures, which resembles the 6c-plane spin spiral in TbMnOa (g5=0.297r) [5^, "6]. 

We trace dynamics of the Mn spins by numerically solving the Landau-Lifshitz-Gilbert 
equation using the fourth-order Runge-Kutta method [27]. The equation is given by 

f ^...f, (4) 

where aG(=0.05) is the Gilbert-damping coefficient. We derive an effective magnetic field 
Hf^ acting on the spin Si from the Hamiltonian 1-L as Hf^ — —dl-L/dSi. Considering 
the strongly reduced Mn moment revealed in the neutron-scattering experiment |24| . we set 
the norm of the spin vector |S'^| = 1.4. The system used for calculations is 20x20x6 in size 
with the periodic boundary condition. For the magnetoelectric (ME) coupling, we consider 
—E ' pij ^ —{E ' 7Vij)Si • Sj where pij is given in Eq. (2) [El [H]. This couphng effectively 
modulates the in-plane ferromagnetic exchanges from Jab Si • Sj to {Jab — E • 7Vij)Si • Sj. More 
concretely, the applied £^||±a modulates the spin exchanges as shown in Fig. [D^d). Here 
\7Tfj\ is calculated to be 3.5x10"^^ jj^Cm from the lattice parameters [25] and the observed 
P {^5000 fiC/m^) for HoMnOa with an up-up-down-down spin order [2S]- This means that 
Ea^l MV/cm induces the modulation | A J^^l^l^^a^^^- 1=0.022 meV. Note that pij are perfectly 
canceled out in the equilibrium state, and never contribute to the static ferroelectricity. 

The electromagnon spectrum is calculated as a response to the weak 5-function pulse. We 
trace the time evolution of the net polarization P(t) = Yl<^i jy^Pijif) after application of the 
pulse where Pij{t) is given by Eq. (2). Then we perform the Fourier transformation of P{t) 
to obtain Im e{(jj). The calculated spectrum is displayed in Fig. H^e), which has two peaks at 
Cl;=0.94 and 2.1 THz, and reproduces well the experimental spectrum of TbMnOs [T^ . 

Using this model, we theoretically demonstrate switching of the spin chirality by application of 
an intense optical pulse. Note first that there are two kinds of 6c-plane spirals with different spin 
chiralities, i.e. clockwise and counterclockwise turns as one propagates along the +6 direction. 
Their spin chiralities % point in the —a and +a directions a and x||+a) so that they are 



Figure 2. (a) Time evolution of the spin chirality X=(Xa? Xb^ Xc) after application of the pulse 
with £"0=11 MV/cm, which shows the chirality reversal. Inset shows spin states before/after 
applying the pulse, and time profile of the applied pulse Ea{t). (b) Energy diagrams of the ah± 
and hc± states in the presence and absence of £^||±a. (c)[(d)] Direction of the net spin chirality 
X and c-axis components of the DM vectors id^j on the in-plane bonds for the clockwise 
[counterclockwise] a6-plane spin spiral, i.e. the ah- [a6+] state when £^=0. (e) [(f)] Staggered 
modulation of the local spin chiralities Xij =t ^Xi,j under £||+a in the ah- [a6+] state. 

referred to as hc- and 6c+, respectively. The chirality % is calculated as a sum of the local 
contributions Xi,i+x ^ Si x 5^+^ and Xi,i+y = Si x Si+y as % = ^J2i(Xi,i+x + Xi,i+y)/S^. 
Note that Xi,i+z is zero because of antiferromagnetic stacking in z-direction. The a6-plane spiral 
spin configuration is also possible in TbMnOa although it is slightly higher in energy than the 
6c-plane one without external fields. Again there are two kinds of a6-plane spirals, ah- and a6+, 
which have c and x||+c, respectively. 

3. Results and Discussion 

Starting with hc- with a, we apply an intense pulse of electric field £(t) = [£a(t), 0, 0] 
along the a axis. Chirality switching by the application of an intense optical pulse has been 
theoretically proposed in Ref. [27j. It has turned out that the switching processes show highly 
nonlinear behaviors with respect to the strength and shape of the pulse. Here we examine the 
pulse of 

Eait) = £;ocoscjtexp[-^], (5) 

which is different from the pulse with sin ut oscillation in the previous study [27] . Here uj and £"0 
are fixed at 2.1 THz and 11 MV/cm, while the full width of the half maximum of the Gaussian 
envelope is taken to be 0.5 psec [see inset of Fig. [2](a)]. In Fig. [2](a), we display calculated time 
evolution of the a-, h- and c-axis components of %. We find a reversal of the spin chirality x 
from he- (xa<0) to 6c+ (xa>0). 

To understand this phenomenon, let us first consider energies of the four chirality states, i.e. 
ah± and hc±. In Fig. Efb), we summarize energy diagrams in the presence and absence of E. 
When £^=0, the a6+ and ah- states are degenerate, and are higher in energy than the ground- 
state 6c-plane spirals. Application of £||±a lifts this degeneracy, and if it is strong enough, 
either a6+ or ah- becomes the lowest-lying state depending on the sign of E. Under the strong 
£||+a (£||— a), the a6+ {ah-) state becomes the lowest in energy. In contrast, the energies of 
hc± are not affected by E. 

More concretely, when £=0, the spins in both a6+ and ah- states are rotating with 
nearly uniform turn angles, and all of the local spin chiralities Xij—Si x Sj have nearly 
the same amplitudes along +c and — c directions, respectively as indicated by green arrows 
in Figs. El^c) and (d). The local chiralities Xij couple with the c-axis components of the 




Figure 3. (a) Schematic figure for time evolution of the energy structure in the 9 space for the 
chirahty-reversal process, (b) That for the chirahty-flop process, (c) Calculated energy structure 
in the 9 space at £^=0. Inset: Spin chirality % and the angle 9 in the cycloidal spin state. 

DM vectors, =td^j, on the in-plane Mn(z)-0-Mn(j) bonds, whose magnitudes are all equal to 
^ab- Since the signs of id^j are staggered as shown by blue arrows, the total DM energy, 

+ d i^i-\-y • Xi,i+y)i cancels out for both a6+ and ab- states, resulting in 
their degeneracy. When E is applied, the modulated ferromagnetic exchanges Jab — ^Jab and 
Jab + ^Jab cause slight spin rotations as indicated by thin red arrows in Figs. [2](e) and (f) to 
decrease and increase of the spin turn angles. Consequently the local spin chiralities become 
subject to a staggered modulation as Xi.j ^ ^Xi.j- The way of the modulation is different 
between a6+ and a6_, namely depending on the sign of %. When £^||+a, the +Ax^j and 
— Ax^j are arranged antiparallel (parallel) to in the a6+ {ah-) state, resulting in the 

decrease (increase) of DM energy. When this DM energy is large enough under the sufficiently 
strong £^||+a, the a6+ state becomes the lowest in energy as shown in Fig. [2](b). In contrast, 
under the strong E\\—a^ the ah- state becomes the lowest in energy. We note that energies of 
the hc± states never change even under E because the a-axis components of the DM vectors are 
alternately stacked, by which the DM energy cancels out despite the staggered modulation of 
the local chiralities, Xi.j ^ ^Xi,j {\\^)- 

Two key ingredients of this chirality switching are dynamical deformation of potential 
structure under the oscillating E via the above-discussed ME coupling and an inertial motion of 
the spin chirality. In Fig. El(a), we depict a cartoon for time evolution of the potential structure 
in the 9 space for the chirality-reversal process from 0=180° (6c_) to 0=0° (6c+). Here 9 is the 
angle between the chirality % and the a axis as defined in the inset of Fig. [3](c). When Ea{t)>0 
as at tr^l psec, 0=90° (a6+) becomes a new energy minimum, so that the chirality % starts 
rotating towards this minimum. Importantly the chirality does not stop its rotation at 0=90° 
immediately, but passes through that minimum or oscillates around it because of the inertial 
force. The DM interaction and the single-ion anisotropy originating from the spin-orbit coupling 
make the rotation of the spin-spiral plane massive, resulting in its inertial motion. When Ea{t) 
becomes again at t >1.5 psec, the system falls into the minimum at 0=0° {hc-^) if the system 
is within the domain of metastability of 0=0°. The system can also be trapped in the local 
minimum at 0=90° (a6+), resulting in the 90° flop of the spin chirality [see Fig. [3](b)]. 

Here the equation of motion for the chirality oscillation around an energy minimum at 0=0o 
can be written by introducing the effective mass fh as m0 + fhj9 = — where 7 is the 

damping rate determined by aQ. Note that the dimension of rh is not M but ML^ since is 
a dimensionless variable in contrast to a variable appeared in the usual equation of motion. 
We calculate the dependence of energy by rotating the spin spiral plane around the h axis, 
which corresponds to V{9). We display calculated V{9) for E^O in Fig. [3](c). Since the typical 



frequency of the oscillation is a few terahertz, the order of fh is estimated as 10~^^ sec^eV. 

This optical chirality switching requires a strong E (|£'o|>ll MV/cm) because it is necessary 
to lower the a6-spiral state energetically relative to the 6c-spiral one. This threshold value can be 
reduced if we properly choose the materials: For example, the solid solutions Tbi-a^Gd^^MnOs 
locating near the boundary between the ah and he spiral phases [28] are promising candidates. 

4. Summary 

In summary, we have theoretically demonstrated the optical switching of spin chirality by 
intensely exciting the electromagnons in multiferroic TbMnOa. We have revealed that the 
oscillating electric-field component of light activates collective rotations of the spin-spiral plane 
with terahertz frequencies via the ME coupling, and their inertial motions result in the chirality 
reversal. The collective spin motion is crucial to this phenomenon, which supplies energy to 
cross the potential barriers. The mechanism here proposed is generic, and can be expected also 
in other spin-spiral-based multiferroics. 
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